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We present a laser cooling scheme for trapped ions and atoms using a combination of laser 
couplings and a magnetic gradient field. In a Schrieffer- Wolff transformed picture, this setup cancels 
the carrier and blue sideband terms completely (up to first order in the Lamb-Dicke parameter), 
resulting in an improved cooling behaviour compared to standard cooling schemes in the Lamb-Dicke 
regime (e.g. sideband cooling) and allowing cooling to the vibrational ground state. A condition for 
optimal cooling rates is presented and the cooling behaviour for different Lamb-Dicke parameters 
and spontaneous decay rates is discussed. Cooling rates of one order of magnitude less than the 
trapping frequency are achieved using the new cooling method. Furthermore the scheme exhibits 
fast rates and low final populations even for significant deviations from the optimal parameters and 
moreover provides good cooling rates also in the multi-particle case. 
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I. INTRODUCTION 

Cooling of atoms/ions is an important precondition for quantum logic and computation [Tj, [2] , wherein an initial 
ground state preparation of the vibrational state is required, for high resolution spectroscopy beyond the Doppler 
limitations and for the analysis of low temperature phenomena (e.g. BEC condensates). Cooling by light (laser 
cooling) for atoms and ions was originally proposed in J3J0]. For trapped ions, sideband cooling [51(5] is the standard 
cooling method, achieved by tuning the laser frequency to the red sideband of the bare atomic frequency, which 
leads to a resonant cooling contribution whereas heating is only induced by off-resonant carrier and blue-sideband 
excitations. To obtain well-resolved sidebands requires the linewidth of the atomic transition to be much smaller 
than the trapping frequency r <C v (strong confinement). Furthermore, taking into account that off- resonant heating 
excitations increase with increasing laser intensity the coupling Rabi-frequency has to be smaller than the trapping 
frequency fi <C v (weak coupling) [2 . Cooling to lower temperatures was reached using the concept of Raman sideband 
cooling [7] , applicable for atoms [5] and ions [H] [TO] , which allows the realisation of very small cooling linewidths by 
Raman two-photon transitions. Further improvements on the way to reach the ground state are achieved with dark 
state concepts based on eliminating off-resonant excitations via quantum interference. In EIT cooling [H] [12] the 
concept of electromagnetically induced transparency [13] is used to eliminate the carrier transition, which gives the 
main heating contribution in the sideband schemes. It uses a three level system coupled by two lasers, wherein one of 
the lasers (the stronger one) creates the absorption spectrum, a Fano like profile with a zero at the carrier transition 
and a peak at the red sideband transition. Therefore stronger couplings can be realised, no strong confinement is 
required due to the three level nature of the system and a large cooling bandwidth (i.e. frequency range for cooling) 
allows simultaneous cooling of several modes. An alternative approach, based on the Stark shift gate [13], is the Stark 
shift cooling [15] , cancelling the carrier transition for a correctly tuned Rabi frequency, which is of the order of the 
trapping frequency. This is achieved in the dressed state picture, i.e. taking the Stark shift of the laser interaction 
(originating from the zeroth order in the Lamb-Dicke parameter contributions of the resonant coupling laser) into 
account. The robust laser cooling scheme [T5] combines the methods of EIT and Stark shift cooling resulting in a 
dark state for the combination of EIT and Stark shift Hamiltonian contribution and leading to an elimination of both 
the carrier transition and the blue sideband heating transition by interference. 

The scheme presented here also eliminates the carrier and blue sideband transitions in lowest order in the Lamb- 
Dicke parameter, using an atomic three level system, a combinations of four laser beams and applying a magnetic 
field gradient. The application of an inhomogeneous magnetic field to trapped ions resulting in a position dependent 
resonance shift and mediating the coupling between internal and external degrees of freedom (e.g. allowing the 
coupling by microwaves with a non- vanishing effective Lamb-Dicke parameter) was discussed in |17l 118] . For our 
method only the proper choice of the gradient strength and the phase relations of the laser beams form the conditions 
for a pure red sideband interaction up to first order in the Lamb-Dicke parameter. 

The paper is organized as follows: In section[TT| the cooling scheme is presented and the influence of the magnetic 
field gradient is analysed. It turns out, that in a Schrieffer- Wolff transformed picture, both the carrier as well as 
the blue sideband transition are eliminated by interference and thus cooling to the vibrational ground state in zeroth 
order is achieved. In section [TTl] the resonance condition for the cooling process is derived. Moreover it is shown that 
the maximal cooling amplitude is limited by the trapping frequency, a fact that arises from the amplitudes of the 
dressed states, which are tuned to resonance. Section |TV] presents an analytic formula for the cooling rate obtained by 
the method of adiabatic elimination and valid in the weak coupling case. Indeed a comparison to numerical results 
shows a good agreement in the range of small Rabi frequencies. In section|V] the dependence of the cooling rate on 
the spontaneous decay rate as well as on the magnitude of the Lamb-Dicke parameter is analysed. Furthermore the 
dependence of the cooling rate and final population on phase fluctuations of the coupling lasers is presented resulting 
in the conclusion that the scheme, albeit becoming slower, still provides fast cooling and low final populations even for 
relatively large deviations from the optimal phase relation. Section |VT| deals with the cooling of multiple vibrational 
modes, i.e. with multiple trapped ions. Only for the center of mass mode a complete elimination of heating terms 
(carrier and blue sideband) is possible; however the strength of the heating terms for other modes is smaller than in 
standard cooling schemes. 

II. COOLING SCHEME AND HAMILTONIAN 

The scheme presented here is based on a complete elimination of the blue sideband heating terms in a Schricffer- 
Wolff transformed picture (see eq. ^ and ^ below) in lowest order in the Lamb-Dicke parameter, which is achieved 
by a magnetic gradient field together with a four laser coupling scheme. It requires an atomic three level system 
consisting of a (fast decaying) excited state \e) and a Zeeman split ground state | + 1) and | — 1) induced by the 
external position dependent magnetic field. Each of the Zeeman split levels is coupled by two lasers with Rabi 
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FIG. 1. Cooling scheme: a) Coupling scheme in the standard basis. Each of the Zeeman split ground state levels is coupled 
to the excited level \e) by two lasers with Rabi frequency Q, Lamb Dicke parameter r\, phase <p and detuning 8. Note that the 
external magnetic field gives, beside the position independent ljo splitting, an additional position dependent energy shift via 
the magnetic gradient contribution, b) Scheme in dressed state basis after Schrieffer- Wolff transformation. See also eq. (fsl . 



frequency Q and opposite Lamb-Dicke parameters r\ [r\ = k y/h/(2 m v) cos 9, with 9 the angle of the laser beam to 
the trap axis and k the absolute value of the wavevector), obtained by irradiating the beams from opposite directions, 
to the excited level |e) having a phase difference ip = n/2 (see figjlj. 

The ion is assumed to be trapped in an harmonic trapping potential with frequency v and for simplicity only the 
motion along the trap axis is taken into account. Moreover, the expansion up to first order in terms of the Lamb-Dicke 
parameter requires the atom/ion to be localized on the scale provided by the laser wavelength (Lamb-Dicke limit, 
rj <C 1). The position dependent magnetic field B(x) has the form 



(1) 

x=0 



and higher terms are not present or can be neglected in the following expansion in terms of the effective Lamb-Dicke 
parameter respectively. The corresponding energy shift due to the magnetic Zeeman splitting is (with the position 

operator x = \J ~^^ ip + b^), wherein and b are the phonon creator and annihilator operators respectively) given 

by 

AE = —jiB = Kuo m a + ?y e // his (6 + & ) m s (2) 

with /i the magnetic moment, m s the magnetic quantum number and Tj e ff = yf=^i7B^ providing a coupling 
between the spatial position and the internal energy, i.e. phonons and internal atomic states. The effective Lamb 
Dicke parameter has a physical meaning as the energy difference between the left and right bound of the ground state 
wavefunction extent in units of hv. The energy shift of the homogeneous part is given by Hujq = —g fis B$, wherein 
g denotes the g-factor and jiB the Bohr magneton. Thus the Hamiltonian of the system (cp. to fig.[l} is given by 
(setting H = 1) 

H =vtf b + u e \e){e\ +Uo (| + 1)(+1| - | - 1)<-1|) + r? e// v (| + 1)(+1| - | - (6 + fit) 

+ Q (|e)(+l| + I + l)(e|) ■ [cos (k x - u L t) + cos (-k x - uj L t + tt/2)} (3) 
+ n (|e)(-l| + | - l)(e|) • [cos (-k x -uj* L t)+ cos {kx-uj* L t + n/2)} 

with u e being the energy of the upper level |e). In order to end up with a time independent Hamiltonian we assume 
for the detuning: 5 = oj e — ojq — u>l = ui e + u>q — ui* L , wherein ujl is the frequency of the first laser pair (coupling | + 1) to 
|e)) and u* L the frequency of the second pair respectively (coupling | — 1) to |e)). In an interaction picture, expanding 
the Hamiltonian up to the first order in the Lamb-Dicke parameter rj and applying the rotating wave approximation, 
the Hamiltonian takes the form 



H =vtfb + S |e><e| + r, eff v (| + 1) (+1| - | - 1><— 1]) (b + ftt) 

^ (1 + i) (|e) (+1| + |e) +r? § (1 + i) (b + fet) (| e > (+1| - | e ) (_1|) + h.c. 



(4) 



In order to cancel the blue sideband heating contributions completely, it is important to adjust the magnetic gradient 
field in a way so that 

V = Veff (5) 
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TABLE I. Magnetic gradients needed to fulfill condition (j5jl for different trapping frequencies. Values are given for the S1/2 to 
P1/2 transitions of Ytterbium-172 ions (369 nm) and Calcium-40 ions (397 nm). In this case the magnetic gradient is given by 
^ = (- 2 g") h fj, B kcos8, wherein 6 denotes the angle between the laser beams and the trap axis (77 = \Jhk 2 /(2 m v) cos 6), which 
in the following is assumed to be zero and g e denotes the electron g-factor. 





172 Yb+ 


40 Ca+ 


Trapping frequency [v/(2 tt)]/Hz 


5- 10 5 


1 • 10 6 


5- 10 5 


1 • 10 6 


Magnetic gradient ff /(T/m) 


607.6 


1215 


564.7 


1130 


Lamb-Dicke parameter r\ 


0.13 


0.09 


0.25 


0.18 



as will be explained later. Values for the necessary magnetic gradient to fulfill condition ([5]) for Ytterbium-172 and 
Calcium-40 ions are given in table|T] 

The basic idea of the cancellation mechanism is that a x (b + b + ) + ia y (b — b + ) generates the red sideband and cancels 
the blue one. Since the phonon coupling introduces a term that is proportional to b + b + and the magnetic gradient 
introduces a term that is proportional to i(b — b + ), we can expect that if the phases are adjusted appropriately we 
can get a complete interference of the blue sideband. 

The effect of the magnetic field is best seen by applying a unitary Schrieffer- Wolff transformation U 19 to the 
Hamiltonian Q 



U = exp{- V (| + 1)(+1|-|-1)(-1|) (6-fot)} 
resulting in the following transformations 

b -+ 6 -„(| + l>< + l|-|-l)<-l|) 
1 + 1) -> \ + l)cxp[- V (b-tf)] 
- l)exp [ 77 (6 — & + )] . 



(0) 



(7) 



Application of the Schrieffer- Wolff transformation to the Hamiltonian Q and keeping only terms up to first order 
in 77 we get (fig.[l]D, note that 77 is equal to t] e ff): 



H' = U HU^ 
Herein we have denned: 



vbH + 5\e){e\ + 



71 



\B) - 



\D) = 



1 

1 

71 



(l + »)|e)<B| + i7nV2(l + *)6|e)(D|+h.c. 

+ 1> + |-1», 
+ 1>-|-1». 



(8) 



(9) 



Note that the first order expansion terms of the Schrieffer- Wolff transformation (%//) of the zeroth order term in 77 
(Lamb-Dicke- Parameter) Hamiltonian contributions (see |4| lead to a cancellation of the first order in 77 blue sideband 
laser interaction contributions under the assumption ([5]). 

The complete time evolution in the Schrieffer- Wolff transformed picture is given by 



at' 



-i[H',p}+£ d p. 



The Liouvillian Superoperator C d describes the spontaneous decay from the excited level \e) to the levels 
I — 1). It exhibits the following form: 

C d = + L d . 

Cq p describes the general spontaneous decay terms in zeroth order in the Lamb-Dicke Parameter: 



E 

i={c+i),(-i)} 



2\i)(e\p\e)(i\ -\e)(e\p-p\e)(e\ 



(10) 
1) and 

(11) 
(12) 
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Cf p is the Schrieffer- Wolff transformation contribution of the spontaneous decay term (in first order in rj) : 

£fp= E 7i(-lYv[(b-^),\i){e\p\em]. (13) 
i={(+i),(-i)} 

Note that recoil only contributes in order r) , leading to a very small correction to the cooling rate (decreasing it) 
due to the dark state nature of the cooling scheme and therefore can be neglected for small r\ (which was tested by 
numerical simulations). This is due to the fact, that the probability for being in the excited state is given by rf 
and the spontaneous decay including recoil from this state is a process in second order in -q and therefore recoil only 
contributes with a strength oc rf . 

If the final phonon population \n — 0) is reached, the system will be decoupled from the laser interaction and will 
end up in the state \D, n = 0) (" dark state"). Note that due to the presence of only red sideband terms the phonon 
state \n = 0) will always be reached (unless ft — 0). In the original (i.e. non Schrieffer- Wolff transformed) picture, 
this final state has the form: 

C/ f \D, n = 0) = \D, n = 0) - X] \B,n = l)+0 {jf) . (14) 

Therefore the presented scheme cools to the vibrational ground state (n=0) in zeroth order in r\. Furthermore the 
Schrieffer- Wolff transformation preserves pure states and thus does not increase the entropy (which remains zero). 
A simple (projective) measurement of the internal state |D) as well as a unitary rotation could further improve the 
achieved final population, that is they would result in a state with n = up to second order in 77. 



III. CONDITIONS FOR OPTIMAL COOLING RATES AND LIMITATIONS 



At first sight it might seem possible to increase the cooling rate by an arbitrary amount just by increasing the laser 
Rabi frequency ft. However a detailed analysis of the ft |e)(i?|-coupling terms (cp. Hamiltonian ([8])) shows that this 
is not the case. These terms are changing the |e)-state (shifting it and mixing it with the |£?)-state) and therefore 
crucially influence the cooling term (~ rjQb\e)(D\) . Optimal rates are achieved for tuning the cooling interaction in 
resonance with one of the dressed states of the ft \ e)(B\ interaction. The Hamiltonian Q in the dressed state basis 
takes the form 

H' =v ¥ b + u u \u) (u\ + \d)(d\+ V nV2(l + i)b {(a%\u)+a d e \d))(D\+h.c.} (15) 

wherein \u) and |d) is the dressed state basis (each having contributions of |e) and \B) states), ui u and ujd are the 
energies of those states (i.e. eigenvalues of the diagonalization procedure) 



w u/d - 2 



i (S± yfS 2 +4 ft 2 ) (16) 



and a" and a e the amplitudes of \u) and \d) in |e) respectively (see also eq. (17)), i.e. the excited level |e) expressed 
in the dressed state basis is given by 

|e) = al \u) + a* \d) . (17) 

In the interaction picture with respect to the first three Hamiltonian contributions, the resonance condition (for 
tuning the \u) state to resonance) is easily seen to be fulfilled by setting lu u — v. This procedure leads to the following 
resonance condition for optimal cooling rates: 

Consequently the detuning is - for a given trapping and laser Rabi frequency - uniquely determined in the optimal 
cooling case. 

To see how the magnitude of the cooling term increases with increasing 51 it is useful to define an "effective Rabi 
frequency" ft e //- If \u) corresponds to the resonant state (the one that has been tuned to resonance), then O e // is 
defined as 

n eff = Cla% (19) 
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FIG. 2. a) \Q e ff\/v is plotted over the Rabi frequency £2 in units of the trapping frequency v. It turns out that the maximal 
achievable |fi e //| is limited to the trapping frequency v (see also eq. ( 20 1 ) . b) Absolute value of |e) amplitudes in the dressed 
states; the red curve represents the resonant state. The minimal amplitude determines the dynamics for a given fi e //. 



and is a direct measure of the strength of the resonant (and therefore relevant) cooling contribution, 
fulfilling condition ( 18 1 it follows that 



v\0\ 



eff\ 



In case of 



(20) 



(see figure Ka). From equation (20) it is clearly seen, that it is impossible to increase the cooling term contribution 



above a limiting upper value, or more precisely, f2 e // is limited to a maximal value, which corresponds to the trapping 
frequency v (\£l e ff\max — v). Therefore the trapping frequency is responsible for an upper limit to the maximal 
achievable cooling rate. 

A second limitation arises from the decay rate 7 and the contributions of the excited state |e) in the new dressed 
state basis. Note that spontaneous decay to both dressed states and the \D) state takes place. The cooling process 
relies on the decay to the \D) state, and thus - to keep the process going - it is necessary that the other states decay 
fast enough, i.e. have a sufficient large |e) contribution. Otherwise the system will be trapped (at least for a certain 
amount of time) in one of the dressed states and thus the cooling rate is reduces significantly. For equal decay rates 
7 in (|T2[>, ( 13 >, equal \e) contributions in the dressed states would lead to the optimal case when taking only the 
spontaneous decay influence into account. The |e) -amplitudes in the dressed states for the resonance case (18 1 are 
plotted in figure^. Note that the minimal amplitude determines the dynamics. After a point of equal (e)-amplitudes 
in the dressed states, the minimal |e) -amplitude decreases again for further increasing f2 e //) which sets the optimal 
rate to a value £l e ff < v i depending on the relation between 7 and v. 



IV. RATE EQUATION FORMULA AND COMPARISON TO THE NUMERICAL RESULTS 

The procedure described in [7J allows to calculate an analytic formula for the cooling rate by adiabatically 
eliminating the internal degrees of freedom under the assumption that the cooling dynamics is much slower than the 
internal dynamics. Moreover it is assumed that the trapping frequency is large compared to the cooling dynamics 
(i.e. r/Q <^ v), which allows a projection on the phonon |n)(n| subspace; i.e., a derivation of a rate equation. This 
second condition is only fulfilled for small values of f2 in our case (otherwise the cooling term is of the order of the 
trapping frequency or for further increasing f2 even bigger). The rate W can be written in the form W — A- — A + 
and the final steady state population is given by (n) ss = A+/(A— — A + ), wherein A + /A—- arc the transition rates 
for increasing or lowering the phonon number respectively. Those quantities consist of two time correlation functions 
of internal operators, which can be calculated using the quantum regression theorem and the optical Bloch equations. 
It follows that, due to the fact that there exist no heating terms in first order in the Lamb-Dicke Parameter, A + is 
equal to zero and thus W = A_ and the final phonon population is given by (n) ss = 0. 
Following the procedure given in [2D], we get for the rate W (see appendix for the calculation) 

W= 8 ^ U2 „ (21) 
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FIG. 3. Cooling rate W as a function of the Rabi frequency Q both in units of the trapping frequency v for r\ — 0.1 and 7 = 5u 
(left plot) / 7 = 0.5 v (right plot). The detuning S for a given Q is chosen to fulfill the resonance condition eq. (18 1. The blu e 
line is the rate obtained by numerical simulations, the red dashed line the one calculated using the rate equation formula ( 23 1 . 



where the rate W is defined via: 



d* W 



-W(n) 



(22) 



This rate (21) gets optimal if the term in brackets is equal to zero and the resulting condition is equivalent to the 



resonance condition given in eq. (18). In that case, the rate takes the following form: 



W 



7 



(23) 



The fact that the rate is proportional to 1/7 is a direct consequence of the adiabatic elimination in high 7. More 
precisely, the rate W is proportional to the spontaneous decay rate 7 of the excited level times the probability p u for 
achieving this level. The latter quantity is proportional to I/7 2 , i.e. the probability for reaching the excited dressed 
state (here: \u)) decreases with increasing dissipation. Therefore the rate is given by: W cx (l/7 2 )7 = 1/7- However 
it does not mean that the optimal achievable rate is getting slower for increasing 7, because it is possible to reach 
higher final values for D, due to the fact, that the (internal) dynamics is less limited by the decreasing |e)-amplitude 
of the resonant state, which occurs for increasing SI beyond a point of equal dressed state amplitudes. Figure[3] shows 
the rate obtained by numerical simulations and that one calculated using eq. (23) for different Rabi frequencies. As 



expected, the rate formula describes the rate quite well in the limit of small 51 values whereas it fails completely in the 
opposite limit (when the cooling contribution is getting stronger and the conditions used for the formula derivation 
are no longer fulfilled). 



V. SIMULATION AND RESULTS 



Figure |4] shows optimal cooling cycles for different spontaneous decay rates (7 = 0.5^ < v and 7 = hv > v). In 
the first case (7 < the rate is of the magnitude of 7) the optimal cooling rate is achieved for £l e ff ~ v (cp. to 
eq. ( [Is] )) , i.e. the case of equal |e) -amplitudes in the dressed states. That is what one would expect, because with 
further increasing fi e // to its maximal value v there is one dressed state with a small |e) amplitude and therefore the 
slower spontaneous decay from this state leads to a decreasing rate despite the increasing strength of the cooling term 
contribution. In the 7 = 5^, 7 > z^, W ^ (cooling rate) case, the spontaneous decay rate is much faster than the 
cooling dynamics. Therefore, as expected, it is possible to increase fi e // beyond the "equal amplitude" situation and 
thus a slightly higher effective Rabi frequency can be realized. The spontaneous decay time of the dressed states is 
much less limiting in this case. Note that this higher is needed to compensate the influence of the quantum Zeno 
effect, which reduces the rate compared to the 7 = 0.5 case for a given f2. The optimal achieved rates W (see also 
eq.p2]» are 0.08 z; {Q = 0.85i/,fi e // = 0.65 v) for the 7 = 0.5 z; case and 0.09 v (fl = 2.8v,Q eff = 0.94 z;) for the 
7 = 5 v case respectively. Rates for different values of Jl are presented in figure [3j 

The dependence of the cooling rate W on the spontaneous decay rate 7 is plotted in fig. [5} As discussed before, 
for higher spontaneous decay rates higher values of Jl are possible because of less restrictions due to the spontaneous 
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FIG. 4. Average population (n) vs. time in units of the trapping frequency v for optimally chosen cooling parameters. The left 
plot describes the optimal case for 7 = 0.5 v (Q = 0.85 v, 8 = 0.28 v, rj — 0.1), the right plot for 7 = 5 v (Q = 2.8 v, 5 = —6.84 v, 
r\ = 0.1). The insets show the final cooling stage of the corresponding figures using a logarithmic scale for the population, 
which clearly reveals the final population of the order of rj 1 . 




FIG. 5. Cooling rate W in units of the trapping frequency v for different values of the spontaneous decay rate 7 (in units 
of the trapping frequency v). The inset shows the optimal value of SI (left) and Q e ff (right) for a given 7. While the rate 
increases fast in regions where the cooling is limited by the spontaneous decay rate (small values of 7) it saturates for high 
spontaneous decay rates, where the rate is limited by the limitation of Q e // (see also fig.[2j| . Values have been calculated for 
77 = 0.1, optimal values for the Rabi frequency and the detuning chosen according to the resonance condition eq. ( 18 1. 



decay and therefore the rate increases with increasing 7. For small values of 7, when the cooling dynamics is limited 
by the spontaneous decay rate and not by the fl e ff coupling, the rate increases significantly with increasing 7. For 
high values of 7 the rate converges to a limiting value for the rate. In this case the cooling dynamics is clearly 
limited by the fact, that f2 e // reaches its maximal value v (the trapping frequency (20), see right inset of fig. [5] for 
the corresponding fi e // values). Thus further increasing Jl in the region Sl e // < v for increasing 7 (which is indeed 
necessary to compensate for the higher quantum Zeno effect) leads to a very small increasing of £l e ff and therefore 
leads to a saturation of the rate with f2 e // converging to its maximal value. 

The dependence on rj (for otherwise fixed parameters) is shown in fig. [6] For 7] <C 1, i.e. in the regime where 
heating terms of higher orders in r] can be neglected, the rate scales approximately linearly with increasing Lamb- 
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FIG. 6. Cooling rate W in units of the trapping frequency v for different values of the Lamb-Dicke parameter r\. The 
corresponding dependence (for optimal rates) is shown in the inset. For small values of r\ (when the influence of higher order 
heating terms is small) the rate scales approximately linear in r\. All values have been calculated for a fixed spontaneous decay 
rate 7 — 5u and optimal detunings according to the resonance condition ( 18 1, including terms up to second order in 77. 



Dicke parameter 77. Note that the rate equation ([23]) includes a quadratic contribution of 77, but also the optimal Rabi 
frequency depends on il, scaling approximately oc lf-y/fj and thus leading to an approximately linear scaling of the 
rate dependence on 77. Thus decreasing 77 leads to increasing il for the optimal condition, which can be explained as 
compensating for the quantum Zeno effect although it deteriorates the spontaneous decay amplitude behaviour. By 
increasing 77 to higher values it reaches some optimal point and decreases again afterwards due to the fact that the 
heating terms in second order (and by further increasing also in higher orders) in 77 become more and more important. 

The cooling scheme concept is designed for the Lamb-Dicke regime (77 <C 1), because only in this regime higher order 
contributions can be neglected. Numerical simulations for higher Lamb-Dicke parameters show that the (minimal) 
final population always increases with increasing Lamb-Dicke parame ter. While the final population for 77 = 0.1 
is given by (n) = 0.01 corresponding to rj 2 as expected according to (14), it is given by (n) ~ 0.06 for 77 = 0.2, 
(71) ~ 0.64 for 77 = 0.3 and (77) ~ 4.5 for 77 = 0.6 (f2 = 13 v, 7 = 5z/). Therefore only in the Lamb-Dicke regime the 
goal to achieve the vibrational ground state can be achieved. With respect to the new final population the cooling 
exhibits for higher Lamb-Dicke parameters in average again an exponential decay, leading e.g. for the above described 
rj = 0.6 case to a cooling rate of W ~ 0.006 v. Note however that the cooling rate and the final temperature are 
not optimized simultaneously for the same parameters, which especially for high Lamb-Dicke parameters leads to 
significant differences in the final populations and rates. The fact that both the final population and rate are worse 
for increasing Lamb-Dicke parameters underlines the restriction to the Lamb-Dicke regime. 

Robustness to phase fluctuations — The presented scheme relies on destructive interference and thus may be fragile 
to dephasing, i.e., fluctuations in the phase difference between the laser pairs around the optimal 7t/2 value. The 
numerical analysis of the phase fluctuation is shown in fig. [7] It can be seen that even at fluctuations which are larger 
than 20% the population is still less than 10 -2 and the rates decrease only by a factor of 2. 

The proposed cooling scheme is faster than previously proposed continuous cooling schemes. The fastest continuous 
scheme that was implemented is the EIT cooling scheme [12] which cooled at a rate of 10~ 4 ta The robust cooling 
scheme |16j can in principle reach a rate of 10 _2 f. In contrary the proposed cooling scheme can reach rates which are 
one order of magnitude less than the trap frequency. 



VI. MULTIMODE COOLING 

The presented cooling scheme was tested for the case of several trapped ions. Furthermore the cooling scheme, 
i.e. the laser coupling of internal and external degrees of freedom as presented in sectionlTT] is applied to one of the 



More exactly, numerical simulations show that the approximate scaling is given by W oc rj 28 and fi oc rj 
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FIG. 7. The left plot shows the final population (n) ss , the right plot the rate W in units of v for different phase differences ip 
between the two coupling laser pairs. Note that ip = n/2 corresponds to the optimal case allowing the complete elimination of 
heating terms in first order in r\. The simulation is performed in the Schrieffer-Wolff transformed picture and an Hamiltonian 
expansion up to second ord er in the Lamb-Dicke parameter with the parameters Q = 2.8 u, 5 — —6.84^ (according to the 
resonance condition eq. ( 18 1), r\ = 0.1 and 7 = 5 v is used. Note that the final population corresponds to the population in the 
transformed picture. 



trapped ions (named ion 1 in the following). It is advantageous to express the Hamiltonian in terms of normal modes, 
wherein the creator and annihilator of the n-th normal mode is denoted by 4 and a n respectively. In contrast, the 
ones of the i-th local mode are denoted by bi and b\ . Expanding equation Q to the multi particle case leads to the 
Hamiltonian 



H=6\e){e\+J2^al l a n +r,u [\ + - | - 1)(-1|] (&i + 

n 

(1 + i) \e){B\ + r} SL(i + i ) \e)(D\ (b, + b{) + h.c. 



(24) 



-^(l + ,)|e>( fl | + ^. 

Herein v n describes the frequencies of the corresponding normal modes, the summation over n equals a summation 
over all normal modes and all internal states correspond to those of ion 1. With the transformation Xi = J2 n Qn 
between the local coordinates and the normal coordinates {Q n } the local operators in (24) can be expressed in 
terms of normal ones using the relation (xj = \J hj {2mv) (pi + b\)) 



(25) 



and thus the Hamiltonian ( 24 ) reads 



H =5 |e) (e| + £ [v n 4 a n + V ' n v [| + - | - 1) (-1|] (a n + 4)] 



^(l + i)\e)(B\ + J2v' n ^Ci- + i)\e}(D\ (a n + 4) + h.c. 



(26) 



wherein rj' n —rj M™ \fv/v n . 

Again, the effect of the magnetic gradient terms is best seen by moving to a Schrieffer-Wolff transformed picture in 
which the magnetic gradient displacements of the individual normal mode oscillators are removed. This transformation 
is given by 



U' = exp 



-X>n[| + l)<+l|-|-l)<-l|](a n -4) 



= J] exp [-fjn [| + 1)(+1| - I - 1>(-1|] {a n - 4)] 



(27) 
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FIG. 8. Cooling of three trapped ions. The line colours blue, red and green correspond to the single modes with increasing 
frequency in this order. The simulation was performed in the original (non transformed) picture using a second order in r\ 
Hamiltonian expansion and the values Q = 2.8^,7 = 5i> and r\ — 0.1. The insets show the corresponding cooling processes 
with a logarithmic scaling of the population. The left plot shows the cooling of the three modes all starting with an initial 
population (n) = 2. The detuning 8 was chosen in a way to set the red sideband terms of the first and second mode to resonance 
respectively, corresponding to continuous and dashed lines. Continuous lines describe the situation when the laser is tuned 
to resonance with the red sideband of the center of mass mode and dashed lines the analogous situation with the breathing 
mode tuned to resonance (see resonance condition eq. (18l). The right plot shows the heating of the excited modes during the 
cooling process of the center of mass mode, which is set to resonance. The final populations are: (ni) = 0.005, (112) = 0.031 
and <n 3 ) = 0.14. 



with fj n = T]' n vjv n — 77 M" {vjv n )' i l' 1 . Note that this transformation differs from the transformation used in sect.[ll| 
what can easily be checked by inserting (25) into ([6]). Applying this transformation (cp. to eq. Q) to the Hamilto- 
nian (p4] leads to H' = U' HU* 



n 



H' =8 |e)(e| + 2J v n a n a n 

+ ^(l + 0|e>ME^V^ 



{l+i)\e)(B\ 





V 




V 


aCj +h.c. 


( 


1 + — 


a„ + 


1 




v n _ 




Vn_ 



(28) 



Consequently, the amplitudes of the blue and red sideband terms for a specific mode n depend on the relation of the 
trapping frequency v to the normal mode frequency v n . A complete cancellation of blue sideband heating terms is only 
possible for the center of mass mode (COM), because only in this case vcom — v - For higher modes heating terms 
are always present, although they are off-resonant by tuning the laser frequency to the red sideband. In general, for 
increasing v n the heating amplitude increases and the cooling amplitude decreases, which leads to equal amplitudes 
for v n — > 00. The latter case corresponds to the situation in cooling schemes that do not cancel blue sideband 
terms, e.g. EIT cooling and sideband cooling. Therefore the presented cooling scheme results in an improved cooling 
behaviour with respect to those schemes also in the multiparticle case albeit not being able to cool other modes than 
the COM-mode without heating terms. 

The cooling of the individual modes (COM mode (blue), first excited mode (breathing mode, red) and second 
excited mode (egyptian mode, green)) for the case of three ions is presented in fig.[8j For three ions 2 M 1 = 
(1/V3, — l/V^i 1/V6)j v x — v> v i — y/3v, V3 — y29/5 v. The left plot shows the cooling starting with an initial 
population (n) — 2 for each mode. Due to the fact that in the higher excited modes heating terms are always present 
while in the COM mode those contributions are eliminated completely in first order in 77, the COM mode is always 
cooled with a faster rate also in the case when the cooling terms of any of the excited modes is set to resonance, i.e. 
the COM cooling transitions are off-resonant. The right plot shows the heating of the non-resonant excited modes 
originally in (n) = during the resonant cooling process of the center of mass mode starting from (ni) = 2. The 
modes are cooled/heated to final populations of (ni) = 0.005, (n 2 ) = 0.031 and (n 3 ) = 0.14. That is the heating 



turns out to be very small and increases with v n as expected from eq. ( 28 ) 
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VII. CONCLUSION AND COMMENTS 



We presented a cooling scheme in the Lamb-Dicke regime, which completely eliminates the influence of the blue 
sideband heating transitions as well as the carrier transitions in lowest order in the Lamb-Dicke parameter, resulting 
in a good cooling behaviour and providing cooling rates of the order of one magnitude below the trapping frequency. 
Albeit the Lamb-Dicke regime required for the Hamiltonian expansion and the resonance condition limit the cooling 
rate well below the trapping frequency, the achieved rates are faster than in comparable existing schemes working in the 
same regime (e.g. EIT cooling [TT], Robust cooling scheme [15]). This is achieved by a proper choice of the magnetic 
gradient strength whose influence was analysed in a Schrieffer- Wolff transformed picture. Optimal cooling conditions 
were derived together with a formula valid in the weak laser coupling case. The dependence on the Lamb-Dicke 
parameter as well as on the spontaneous decay rate was presented. Moreover the cooling behaviour under deviations 
around the optimal phase differences between the coupling laser pairs was analysed showing that deviations of about 
20 percent only reduce the rates by a factor of two and the final population is still of the order of i] 2 , both very 
satisfying values compared to other cooling schemes working in the Lamb-Dicke regime. Furthermore an analysis of 
the multiparticle case revealed good cooling behaviour also for several modes albeit not being able to eliminate the 
carrier and blue sideband terms completely up to first order in 77 in other modes than the center of mass mode. 
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Appendix A: Rate equation derivation 



In this appendix, the calculation leading to the rate equation formula (21) is presented. The master equation (10) 
up to first order in 77 can be written in the following form 

dp 



dt 



= C aP +(C\+C\) p 



(Al) 



wherein the lower index indicates the order in the Lamb-Dicke parameter. The Liouvillians are given by 

Cq p = (£ ,ext + £a,int) p= -i [H ,ext, p] + [-i [H 0>i nuP] + £q p] ( A2 ) 



C{p = -i[H llP ], C\p = Cip 



(A3) 



with the definitions of the spontaneous decay contributions Cq p and Cf p given in eq. (12) and eq. (13) respectively. 
The complete first order Liouvillian C\ is defined as the sum of the two first order contributions, i.e. C\ = C\ + C\. 
The zero and first order Hamiltonian contributions (cp. to eq. ([8])) are 

H ,ext = vtfb 

(l + i)|e)<B|+h.c 



ffo,mi = 61 |e)(e| + 



V2 



(A4) 
(A5) 



Hi=br)nV2(l + i)\e){D\+h.c.= ^ R n F n 

"£{1,2} 

where R n are pure external operators and F n pure internal ones: Rl :=b, R 2 :=R\, Fi := T) ft y/2 (1 + i) |e) (D\ and 

f 2 = fI 

The projector V on the relevant subspace is given by 



Vp = V e V i p = np ss 
V i p= lime CoMt p = psstnip) 

t—¥oo 
00 

V e P = Y. \n){n\(n\p\n) 

n=0 



(A6) 
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with the external/ vibrational density matrix contribution p = V e tvi{p), the internal steady state p ss — \D)(D\ and 
V 1 and V e the projectors on the internal steady state (assuming that the internal dynamics is much faster than the 
cooling dynamics) and the vibrational \n)(n\ subspace respectively, trj denotes the trace over the internal part. The 
master equation for the relevant external part V p up to first order in the Lamb-Dicke parameter is given by (see e.g. 



) 



_d 
df 



poo 

Vp = Vd (-Cu l )C 1 T P = / m>£xe Cot C 

Jo 



Vp. 



Tracing equation (A7| over the internal part, the master equation for the external part p is obtained 
d 



df 



dt V e tT! (£ ie £ot £iPp 



&t V e tn (C\e Cot C\Pp 



(A7) 



(A8) 



wherein it was used that tr i(Vp) = V e trj(p), which follows directly from eq. (A6). The last step takes into account 
that the L\ terms (originating from the Schrieffcr- Wolff transformation of the spontaneous decay terms) do not 
contribute in the master equation, because 



^p Pss = 0, Z x p ext pi 







with pss = \D ) (D\ and furth ermore assuming equal 7-rates in L\. 
Inserting (A5| into eq. (A3) the first order Liouvillian C\ can be written as 



using the definitions £™ ext p = [Rn, p] and C" int p = [F n ,p]. Thus, inserting ( A10) into ( A8) and noting that 



one obtains 
d 



tT l(£lint Pext <S> Pint) = Pext ^ l{[F n , Pint]) = 



PPs 



(A9) 



(A10) 



(All) 



f-OO 

■■(-i)J2 atV e tv I {F n Cl ext e c ° t C\^ 

n J 

= (-l)E / dt 'P eC le X te Co ^ t {ci ext ptT I (F n e Co ^ t F n , P ss)+pRn>tr I (f„ e £ °— *£^ nt p ss ) } (Al2 ) 

n.n' 

= ("I) E r d ^ e£ I^ e£ °' ea=t ' { C lext P (Fn(t) F n ,(0))ss + IX Rn- ([F n (t), F n , (0)])„} , 
n.n' 

with the internal operator two time-correlation functions (F n (t) F n > (0)) ss = tr/ (F n e CoAntt F n i p ss ) and (F n / (0) F n (t)) ss 
tr/ (F n e c ^ 1 Pss F n ,). 

Furthermore, using that [exp(C . ex t t),p] =0 and defining R n (t) = e £o <!Xt 1 R n we end up with 



d f°° 

-,i= (- 1 )Y,J dtVC {[Rn,[RrAt),p]} (F n (t)F n 



(0))« + [Rn, H Rn' (*)] [F n (t], F n , (0)] } 



(A13) 



Inserting the concrete expressions for the external and internal operators (see definitions below eq. (A5|) and applying 
the projection V e (eq. (A6|) the master equation for the external part turns out to be 



dt 



wherein 



t-i = S 12 (-v) (6 f pb-btf p) + S 2 \[y) (bptf -bHp) +h, 



Sii{v)= / dte^(F(t)F 3 (0)) s 



(A14) 



(A15) 



The time evolution of the vibrational state population {p — X]^Lo-Pn |^)("|) following from eq. (A14| takes the form 
^-p n = (n\^-p\n) = (n+ 1) A-p n+ i - [{n + 1) A + + nAJ\ p n + nA + p n ^i (A16) 

at at 
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with the phonon number increasing / decreasing rates 

A- =2Re(SaiM), A + = 2Re(Si 2 (-t/)) . (A17) 
Therefore the equation of motion for the phonon population expectation value (n) is found to be 



d d d 

dt ^ = dt ^ Pn n = ^ ™ di P " = ~ A+) ^ + A+ (A18) 

n=0 n=0 



with the steady state solution (n) ss 



A_-A, 



(A19) 



Thus it remains calculating the rates A- and A + , i.e. calculating the two-time correlation functions Sij(v) ( A15 1 . 
Using the Bloch equations (the internal equation of motion in the zeroth order in the Lamb-Dicke parameter) and 
the quantum regression theorem, one obtains the results in section |TV) 
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